Scale-free patterns at a saddle-node bifurcation in a stochastic system 
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We demonstrate that scale-free patterns are observed in a spatially extended stochastic system 
whose deterministic part undergoes a saddle-node bifurcation. Remarkably, the scale- free patterns 
appear only at a particular time in relaxation processes from a spatially homogeneous initial condi- 
tion. We characterize the scale-free nature in terms of the spatial configuration of the exiting time 
from a marginal saddle where the pair annihilation of a saddle and a node occurs at the bifurca- 
tion point. Critical exponents associated with the scale- free patterns are determined by numerical 
experiments. 

PACS numbers: 82.40.Ck, 82.40.Bj, 05.40.-a, 05.65.-|-b 



Scale-free patterns are widely observed in nature. The 
most familiar one in statistical mechanics might be that 
in particle configurations at the gas-liquid critical point 
The patterns are characterized by a correlation 
length which is given by a power-law function of the dis- 
tance from the critical temperature. Another type of 
scale-free patterns was found in diffusion-limited aggre- 
gation [^l and spinodal decomposition [3|. In this type, 
the length scale that characterizes the patterns grows as 
a power-law function of time, without the fine-tuning of 
the system parameters. 

Recently, a new type of scale-free pattern called dy- 
namical heterogeneity has been found in glassy systems 
0, i i, S i, i [H • Such patterns are fascinating because 
they become visible only by quantifying some dynamical 
events (called bond-breaking events or unlocking events) 
during a particular time interval . Let us consider dense 
colloidal suspensions as an example. We denote the dis- 
placement of particles in a region around a position x 
during a time interval t as q{x,t). Then, the correla- 
tion length £^{t) of the pattern q{x, t) becomes maximum 
at a time and ^(i*) exhibits divergent behavior as a 
function of the distance from a critical parameter value. 

The characterization in terms of the time-dependent 
correlation length ^(i) is similar to that of growing pat- 
terns, while the existence of the critical parameter value 
is similar to critical phenomena in equilibrium systems. 
Such a coexistence of these features implies that q{x, i*) 
indeed belongs to a new type of scale-free patterns. Since 
its nature is simple and nontrivial, we expect that there 
exists a wide class of systems that exhibit such a type of 
patterns. Motivated by this expectation, we propose the 
simplest model among them. We then intend to elucidate 
the nature of the new type of scale-free patterns. 

Model: In the present paper, we study a spatially ex- 
tended stochastic system that exhibits relaxation behav- 
ior at a saddle-node bifurcation. Although there is a 
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certain idea that connects this model with the under- 
standing of glassy systems [III [3i <io consider 
the relationship in this paper. Instead, we regard this 
model as a typical system that undergoes an elementary 
bifurcation under the influence of noise. By elementary 
bifurcation, we imply pitch-fork bifurcation, Hopf bifur- 
cation, and saddle-node bifurcation [l^. Among them, 
the first two bifurcations under the influence of noise have 
been studied intensively in the the context of critical phe- 
nomena [13, HE 13 ■ Therefore, a saddle- node bifurcation 
with noise might be related to a new class of critical phe- 
nomena. 

Specifically, we study a coupled Langevin equation in 
a one-dimensional lattice = 1, 2, • • • , iV}. Let (pt be 
a one-component quantity defined at the z-th site, where 
we assume periodic boundary conditions (pa ~ (j)N and 
(j)N+i = 4>i- Then, 0^ obeys 

^t(|)^ = /(0,) + k(?!)^+1 + (/).-1 -2(/),)+^«, (1) 

where = — (/)(((^— l)^-|-e) with small e; k is a coupling 
constant and S.i represents Gaussian white noise that sat- 
isfies {it{t)^jit')) = 2T6{t - t')5^j. The noise intensity T 
is assumed to be a small positive constant. The potential 
function f ((/)), which is defined by / = —d,f,v, has a single 
minimum at = when e > 0, while a pair comprising 
the minimum and maximum appears around 0=1 when 
e becomes negative (see Fig.[T]). The minimum and max- 
imum correspond to the node and saddle point, respec- 
tively, in the deterministic equation dt(j) = f{(j))- Such 
a qualitative change in the trajectories at e = in this 
deterministic equation is called saddle-node bifurcation. 
We call the fixed point <f> — 1 for e = marginal saddle 
because this saddle vanishes when e > 0. A characteristic 
feature of ([1]) is that even small noise drastically affects 
the trajectories passing through the marginal saddle. For 
example, when e < 0, the deterministic trajectories start- 
ing from (j) = 1.2 reach the fixed point near <j) — \, while 
trajectories under the influence of noise escape from this 
fixed point and finally arrive at the region near the glob- 
ally minimum point </> = 0. In the argument below, in 
order to extract a simple relation, we focus on the case 
e = at which the saddle-node bifurcation occurs. 
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FIG. 1: Potential w as a function of (j> for several values of e. FIG. 3: Spatial pattern at time t = 828 for the system with 

T = 10"^ N = 4096. 



Since the model we study might be the simplest in the 
class of spatially extended stochastic systems in which 
a saddle-node bifurcation occurs locally in space, our 
model has some relevance to many experimental sys- 
tems. One example is fluctuating motion of a single 
one-dimensional polymer that is subjected to a trapping 
potential in a liquid (see Ref. |i3|)- Another example 
might be found in stochastic reaction-diffusion systems 
because a saddle-node bifurcation is observed in chemi- 
cal reactions (see Ref. [3 ^ ^'^ illuminating example). 
Although multiplicative noise generally occurs in chem- 
ical reaction systems, we expect that the multiplicative 
nature does not affect the behavior reported below unless 
the noise intensity becomes zero at the marginal saddle. 

We investigate ((IJ by numerical simulations. Con- 
cretely, we employ an explicit discrete method with a 
time step At — 0.01, where the error is 0(At^/^). We 
have confirmed that the important results reported be- 
low are kept to be valid when we select Ai — 0.001. We 
also restrict our investigations to the case k — 1 We ob- 
tained the same result for the case k = 0.5. That is, a 
fine-tuning of the parameter value is not necessary [l9|. 
Finally, since we are interested in the patterns emerging 
from a homogeneous state, we assume that 0i(O) = 1.2. 
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FIG. 2; Contour curves defined by (/>i(t) = 0.5 for several 
values of T. iV = 4096. 

Preliminaries: In order to observe the time evolution 
of 0i, we display contour curves defined by (/)i(t) = 0.5 
for several values of T (see Fig. Each contour curve 
distinguishes the late stage ((/)i(<) — 0) from the early 
stage {(j)i{t) ~ 1) because the two stages are connected 



in a short time interval around the exiting time from the 
marginal saddle at each site i. Thus, for example, the 
pattern 0^ at time t = 828 for the system with T = 10~^ 
consists of domains 0^ ~ 1 and (pi ~ 0, as shown in Fig. [31 
where the time t = 828 is selected such that the spatial 
average of (pi becomes 0.5. It should be noted that such 
a pattern appears only around a particular time t — t^. 
Furthermore, the typical length scale of the pattern at 
the time t* increases as T is decreases. This suggests 
that (f>i{t^,) is scale- free in the limit T 0. 

The behavior of (pi{t) is similar to that of a field q{x, t) 
describing a dynamical event in glassy systems (see Ref. 
Q). In order to make the analogy more explicit, we con- 
sider the spatial average of (pi, (p{t) = J2i4'i{t)/^- In 
Fig. m we show the ensemble average {4>{t)) for the sys- 
tem with T = 10~", where a = 3,4, 5, and 6. The staying 
time at the marginal saddle increases as T is decreased. 
The graph of (0(i)) is similar to a time correlation func- 
tion of density fiuctuations in glassy systems. Further- 
more, we measure the fluctuation intensity 

x4t)^N{{p\t))~{Pit)f)- (2) 

As shown in Fig. [5l Xipi^) O'^g peak at a time t — t^, 
which provides a precise deflnition of intuitively used 
in the previous paragraph. Then, Xcl>{t*) increases as T 
is decreased. In fact, the inset of Fig. [5] suggests that 
X0(t*) exhibits the power-law divergence 



(3) 



where the value of the exponent 7^ is close to 1/3. This 
divergent behavior corresponds to the dynamical hetero- 
geneity in glassy systems. These observations naturally 
lead us to a conjecture that the correlation length of 
(pi{t^) exhibits the power- law divergence as a function of 
T. In this paper, we do not investigate 0i(i*) directly but 
propose a novel quantity that characterizes the divergent 
behavior. 

Result: Our basic idea for the characterization of the 
patterns is to employ a special solution (p^, (t) of the equa- 
tion dt(p — —d^v{(j)) under the conditions (p{t) — > 1 for 
t -00, p{t) for t ~* 00, and 0(0) = 0.5. We then 
express the time evolution of the patterns as 



pi{t) = p^{t- 



(4) 
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where the exponents 7 and v characterize the divergences 
of the ampUtude of the spectrum x(fc„; T) and its length 
scale, respectively, in the limit T — * 0. The exponents 7 
and V are determined such that x{f^n', T) with T — 10"*^, 
where a = 3,4,5, and 6, collapse into one universal 
curve. However, after some trials, we find that a log- 
arithmic correction appears in the scaling relation. To 
demonstrate it, we express x{kn]T)T as a function of 
A:„T-i/3/(_inT)i/2 in Fig. El which suggests that the 
scaling relation ([7]) with the logarithmic correction ap- 
pears to be valid for the values 



FIG. 4: (A as a function oil. N = 1024. 
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FIG. 5: x<t> 8-8 a function of N — 1024. Inset: x<ii(**) ^ 
a function of T in a log-log scale. The guide line represents 

where 0i represents the exiting time from the marginal 
saddle, which is defined by (t>i{9i) = 0.5. By this defini- 
tion, the pattern of 9i corresponds to the contour curve 
shown in Fig. [51 For a typical trajectory whose weight 
is large, pi{t — 9i) is expected to be smaller than the 
special solution. Here, we conjecture that the fluctua- 
tion intensity of Oi exhibits a power-law divergence of T 
because the divergence is related to the Goldstone mode 
associated with the symmetry breaking of the time trans- 
lational symmetry. We thus characterize the conjectured 
scale- free patterns by the statistical quantities of 9i. 
Concretely, by introducing the Fourier transform of 6i 

as 

1 ^ 

we measure the spectrum 

x(fc„;T) = iv(|0>„)l'), (6) 

where fc„ = 2Trn/N with n = 1,2, ■•■ ,7V/2. Here, ( ) 
represents the ensemble average. In order to extract the 
singular behavior of x(fc„; T) in the limit T — > 0, we first 
assume a scaling relation 

X(fc„;r) = T-'xiknT-"), (7) 
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FIG. 6: x{kn\T)T as a function of KT-^'^ / {-^InTY^'^ for 
several values of T. N = 4096. The statistical error-bar 
is about the symbol size. By comparing the systems with 
N = 512, 1024, 2048, and 4096, it is found that the finite-size 
effects are fairly small. 

The exponent 7 is related to the scale dimension of the 
quantity 9i. In order to confirm this fact explicitly, we 
consider the spatial correlation function 

[OAn) = r-(^-'')^f]e-^"('-")^^x(if«),(io) 

n=l 

where K,, = KT-'' and AisT T-''2'k/N. From pO|. 
we find a scaling form 

ei=T-^'i-'''^/'^Q{T''l), (11) 

where Q{x) is a fiuctuating field whose distribution func- 
tion is independent of T in the limit T ^ 0. By combin- 
ing ([4]) with this scaling form (|11[) , we arrive at our main 
claim that the scale-free patterns are characterized by 

Ut)-4>*{t-T-^/^Q{T^'H)), (12) 

which clearly indicates the scale-free nature at a saddle- 
node bifurcation point in the limit T ^ 0. We remark 
that the statistical quantities of (j)i , including the critical 
exponent 7^ in ([3]), can be calculated theoretically from 
our result We believe that the representation (fT^ 

captures the essence of the scale-free pattern. 
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Concluding remarks: We have proposed a new uni- 
versality class consisting of stochastic systems undergo- 
ing a saddle-node bifurcation. A remarkable feature of 
this class is that the criticality near a saddle-node bifur- 
cation originates from the fluctuations of the time pass- 
ing through a marginal saddle. We expect that there 
is a rich variety of systems belonging to this class, al- 
though a quantitative experimental measurement seems 
to be challenging. An important example includes neu- 
ronal avalanches [13], which have recently been studied 
extensively. The analysis of a simple mean-field model 
has revealed the critical nature clearly 2l|. We conjec- 
ture that scale-free patterns similar to those studied in 
this paper will be observed in systems related to neuronal 
avalanches. 

Furthermore, we have pointed out that the cooper- 
ative behavior observed in our model is related to the 
dynamical heterogeneity in glassy systems. We wish to 
emphasize that this similarity is not superficial. In fact, 
recently, we have found that the dynamics of fc-core per- 
colation in a random graph exhibits a jamming transition 
via a saddle-node bifurcation of some order parameter in 
the thermodynamic limit p^ . 

The most important result is the expression ([12]), 
which claims that the criticality of the pattern observed 
at a particular time is characterized by 9i . The approach 
based on this idea was addressed in the theoretical anal- 
ysis of a model [1^ . We expect that this approach would 



elucidate the nature of the criticality of jamming systems. 

Before ending this paper, we present some important 
future problems. The first problem is to understand 
the behavior of d-dimensional systems. Theoretically, we 
should start estimating the upper-critical dimension by 
using physical arguments. The second problem is to un- 
cover possible normal forms in multi-component systems. 
Even if the deterministic part exhibits a saddle-node bi- 
furcation, the effect of noise is not determined uniquely 
in such systems. There might be some systems in which 
the exponents are altered from the result we obtained in 
this paper. Finally, the e-dependence of the length scale 
should be clarified. These will define another exponent, 
which might be more relevant in the context of critical 
phenomena. 

In order to consider these future problems, we need 
to develop the theoretical analysis. We have just finished 
the mean field analysis, which will be published soon [231 . 
According to the result, the field 9 scales as ~ y-i/s 
and there are scaling relations as functions of eT"^^^. 
From the naive dimensional analysis, the length scale of 
the system with e ^ turns out to be proportional to 



-1/2 



which leads to 7 = 1 and v = 1/3. 
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